A Study of Meson Correlators at Finite Temperature * 
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We present results for mesonic propagators in temporal and spatial directions at T below and above the 
deconfining transition in quenched QCD. Anisotropic lattices are used to get enough information in the temporal 
direction. We use the Wilson fermion action for light quarks and Fermilab action for heavy quarks. 
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1. Introduction 

Hadronic properties are expected to change 
with increasing temperature, especially near to 
and above the QCD phase transition. Forthcom- 
ing heavy ion collision experiments will give the 
opportunity to investigate them, and require the- 
oretical studies of signals of the phase transition 
and the QGP formation. Lattice QCD simula- 
tions allow us to perform model independent non- 
perturbative studies. 

So far most calculations of hadronic masses at 
finite temperature have handled only the screen- 
ing mass |l]] but not the pole mass, which has a 
fundamental importance. In addition to the in- 
trinsic ambiguity in defining pole masses at finite 
temperature the main difficulty comes from the 
shortness of the "temporal" (temperature) exten- 
sion, l T = 1/T. 

Here, we attack the problem by using optimised 
hadronic operators on anisotropic lattices. Our 
strategy is as follows: 

(1) to use anisotropic lattices to obtain a good 
resolution in the temporal direction. 



(2) to seek for a good operator which has a large 
overlap with the state of interest at zero temper- 
ature. 

(3) to observe temperature dependent effects by 
using this optimised operator on finite tempera- 
ture lattices. 

In the following, we report results on light and 
heavy mesons. Since part of the former have al- 
ready been reported Q2J , we briefly summarise the 
current results. On the latter, we describe our 
computation and preliminary results. 

2. Results on light mesons at finite tem- 
perature 

We use lattices of 12 3 x N T with N T = 72, 20, 16 
and 12 at /3 = 5.68, jg = 4 , in the quenched ap- 
proximation 0. Calibration using Wilson loops at 
N t = 72 gives the anisotropy £ = a a /a T = 5.3(1) 



and 



0.85(3) GeV . Our lattices correspond 



to temperatures T ~ 0,0.93T C , 1.15T C and 1.5T C 
respectively. 

For light quarks, anisotropic Wilson fermions 
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are employed. After Coulomb gauge fixing, the 
PS meson is measured by the following operators; 

G PS {r,t) = w(mi)w(m 2 ) 

< Tr[M(rn 1 ,0;n,t)Af(m 2 ,0;n + r,t) t > (1) 

where M(fn, mo; n, no) is a quark propagator and 
w(m) is exp(—a\rh\ p ), which we call "exp-source" 
here. After tuning the parameters a and p, the 
correlator at T = is dominated by the lowest 
state 0. Then the correlators at finite temper- 
ature are measured: It turns out that the effec- 
tive masses change significantly above T c , while 
T ~ 0.93T C (N t = 20) reveals almost the same 
features as T = 0. 

We also examine the spatial structure of the 
correlator above T c . In Fig. |l| we compare 
the Ps "wave functions" normalised at x = 0, 
Gps(x,t)/G Ps (0,t), at several t for T ~ 0.93T C 
and T~ 1.5T C . Also shown is the free quark case 
(m q a a = 0.1, 7f = £ = 5.3). The exp-exp source 
appears somewhat too broad at T ~ 0.93T C : the 
quarks come close to each other while propagat- 
ing in t. An unexpected fact is that this is also 
the case at T ~ 1.5T C : the spatial distribution 
shrinks and stabilises, indicating that even at this 
high temperature there is tendency for quark and 
anti-quark to stay together. This is in clear con- 
trast to the free quark case which never shows 
such a behaviour regardless of the source. The 
same feature holds for the other mesons at all 
quark masses. 

Fig. shows the temperature dependence of 
"masses" which are extracted from the "exp-exp" 
correlators at the three points [N t /2 — 2,N t /2]. 
Though the values above T c have large system- 
atic errors, pole masses indicate chiral restoration 
there. We also show in the figure the screening 
masses: their difference with the pole masses is 
manifest. It is also noted that mesons degenerate 
approximately irrespective of their spin. 

3. Results on heavy quark systems at finite 
temperature 

Our target is mainly charmonium physics, e.g. 
investigations of the mass shift near T c [|[ and 
J/ip suppression in the plasma phase In or- 
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Figure 1. Ps wave function (n a = 0.086, exp-exp 
source) normalised at r = 0. Also plotted is the 
initial distribution of separations as given by the 
source. T ~ 0.93T C (left) and T ~ 1.5T C (right). 

der to treat heavy quarks with the currently avail- 
able computer power, we use the effective theory 
approach. We adopt Fermilab action with 0(a) 
improvement ||. On the anisotropic lattice, the 
heavy quark operator is 

K(x,y) = S x ^ y - k t {(1-74)7+4 + (1+74)T_ 4 } 
-«a E « r - 7i)r-H + (r + li)T-i] 

i 

-r ^KtCEg<T4iF 4 i(x) + K 8 CB^gcrijFij(x)j 8 X:V . (2) 

where, T ±M = U± f _ l (x)S x± p,. y , n s = l/2(m + £ + 
3r), Kt = £,k s , and r = l/£. Mean- field improve- 
ment is done for ce and cb- 

Introducing the anisotropy results in changing 
the free quark dispersion relation. In eq.(||), one 
notices that a larger anisotropy £ causes a smaller 
Wilson term. Then the problem is how the con- 
tribution from the momentum region at the edge 
of the Brillouin zone becomes larger. In the 
heavy quarkonium, typical energy and momen- 
tum scales exchanged in the meson are of order 
mv 2 and mv respectively. For the charmonium, 
v 2 ~ 0.3, and the typical energy scale is ~ 500 
MeV. Let us consider our lattice: a^ 1 — 0.85 GeV 
and £ = 5.3, m ~ 0.3 roughly corresponds to the 
charm quark mass. From the dispersion relation, 

3 This form differs from that of ]6 



0.4 



0.3 



0.2 



0.1 



0.0 





* o m Ps 




• 




□ m s 


1 


m A 






I 




m 





0.00 0.02 0.04 0.06 0.08 0.10 0.12 
1/NJ 

Figure 2. Temperature dependence of the mass 
7ji (open symbols) and screening mass (full 
symbols) in units a^Mn the chiral limit. The ver- 
tical gray lines indicate T c . The data correspond 
to T ~ 0, 0.93T C , 1.15T C and 1.5T C . 

E(p z = a/ir) - E(0) ~ 0.5 GeV. Though this re- 
sult seems not sufficient, the effect is expected to 
be small for the ground state. With increasing 
lattice cutoff, this problem becomes milder. 

We used the same configurations as for the 
light meson study. The anisotropy £ is fixed 
by using the dispersion relation. At present one 
set of parameters is available: K a — 0.0984 and 
7i? = 3.67 which corresponds to the vector meson 
mass TTiy ~ 3.f0 GeV. 

Fig.|| shows the effective masses of Ps and V 
mesons at finite temperature. As for light mesons, 
on the N T = 20 lattice the lowest vector and pseu- 
doscalar states keep almost the same masses as at 
T = 0. Above T c we observe, as in the light meson 
case, that the V meson effective mass goes below 
the Ps meson. However, the Ps effective mass 
stays at almost the same position as at T = 
while the V effective mass decreases fast, in con- 
trast to the light meson case where both of them 
increase with T. 

4. Summary and Outlook 

For a light quark system, the wave function in- 
dicates a strong correlation between quark and 
antiquark even above T c . The pole masses show 
chiral restoration above T c and a behaviour quite 
different from that of the screening masses. For 
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Figure 3. Effective mass with k„ ~ 0.0984 and 
7f = 3.67. At N T — 72 Pseudoscalar and Vector 
meson mass are 0.6756(7) and 0.6890(8) respec- 
tively. Full symbols show effective masses con- 
structed by free quarks 

heavy quark systems we use the Fermilab action 
on an anisotropic lattice. The effective mass of 
the vector meson drops sharply above T c . Further 
investigations with finer lattices are undergoing. 

The calculations have been done Intel Paragon 
XP/S at INSAM, Hiroshima Univ. This work 
is supported by the Grant-in- Aide for Scientific 
Research by Monbusho, Japan (No 10640272, No 
11440080). 
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